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SPECTRA  FROM  MOLECULAR  DYNAMICS1 
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Following  the  work  of  Roy  Gordon (1)  we  are  accustomed  to  deducing  molecular 
dynamics  from  spectra  Here  we  take  the  inverse  approach  and  compute  spectra  from 
molecular  dynamics.  Specifically,  from  atomic  trajectories  we  compute  infrared,  elec¬ 
tronic,  and  nonresonance  Raman  spectra.  Fig.  1  illustrates  the  general  technique 
From  an  initial  set  of  atomic  positions  and  velocities  and  the  forces  among  the  atoms 
we  calculate  from  Newton's  second  law  the  atomic  trajectories.  For  infrared  or  elec¬ 
tronic  absorption  spectra  we  connect  to  the  radiation  field  through  the  dipole  moment 
ft.  and  for  Raman  spectra  we  connect  through  the  polarizability  tensor  P.  The  tech¬ 
nique  can  be  viewed  as  an  expression  of  classical  linear  response  theory  or  as  an  exer¬ 
cise  in  classical  electromagnetic  theory.  We  use  a  power  spectral  approach  symbolized 
as  D[],  taking  advantage  of  fast  Fourier  analysis  with  windowing  and  associated  correc- 
tions(2, 3).  An  average,  symbolized  by  <>,  of  the  spectrum  is  taken  over  the  ensem¬ 
ble  appropriate  to  the  experimental  conditions  (for  example  a  constant  temperature  and 
density)  and  simple  quantum  corrections  to  the  spectra  are  applied  where  needed 

Molecular  absorption  spectra  can  be  rotational,  vibrational,  and  electronic.  All 
normally  involve  interaction  of  nuclear  position  and  radiation  field  through  the  dipole 
moment.  Fig.  2  shows  an  example  of  a  gas  phase  vibrational-rotational  infrared  spec¬ 
trum.  Elsewhere(2)  we  have  also  computed  pure  rotational  absorption  spectra  and  have 
extended  vibrational  and  rotational  infrared  spectral  calculations  to  the  liquid  phase. 


1  We  thank  the  National  Science  Foundation,  Chemistry,  the  Office  of  Naval  Research, 
Chemistry,  the  National  Aeronautics  and  Space  Administration,  Ames  Research  Center,  and 
the  National  Institutes  of  Health,  Division  of  Research  Resources  for  providing  the  support 
which  has  made  this  work  possible. 
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FIGURE  I.  Calculation  of  absorption  and  scattering  spectra  from  molecular  dynamics. 
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FIGURE  2.  Band  contours  for  the  infrared  gas  phase  CO  fundamental  vibrational- 
rotational  band! 2).  The  triangles  show  our  essentially  classical  molecular  dy  namics  calcula¬ 
tion,  the  solid  line  shows  an  accurate  quantum  calculation  and  the  dots  show  the  experimental 
measurement (4). 
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Fig.  3  shows  the  computed  equilibrium  and  picosecond  transient  electronic 
absorption  spectra  of  I2  photodissociation  in  liquid  Xenon.  The  connection  to  the  radia¬ 
tion  field  is  through  the  transition  dipole  moment,  and  the  Iodine  atoms  are  sufficiently 
massive  that  quantum  corrections  are  negligible.  In  addition  to  one  photon,  or  absorp¬ 
tion  spectra,  we  can  also  compute  two  photon,  or  scattering,  processes,  as  illustrated  in 
Fig.  3  for  a  liquid  solution  Raman  spectrum 
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FIGURE  3.  The  left  panel  shows  the  electronic  absorption  spectrum  for  l2  in  liquid  Xe¬ 
non,  made  up  of  transitions  from  the  ground  stale  to  the  excited  A,  B,  and  B  ''  states! 5).  The 
computed  equilibrium  room  temperature  spectrum  is  shown  at  top,  with  the  gas  phase  l2  band 
contour  as  measured  by  Tellinghuisen(b)  shown  as  dots.  Below  are  the  computed  transient 
spectra  from  0  to  900  picoseconds  for  the  reaction  sequence  beginning  with  photodissociative 
excitation  to  the  A  state,  then  either  /  atom  escape  or  solvent  caging  followed  by  recombina¬ 
tion  to  form  a  highly  vibrationally  excited  /j  molecule  which  loses  its  energy  to  the  solvent. 
The  right  panel  shows  the  isotropic  and  anisotropic  Raman  spectra  of  S’2  in  liquid  Ar.  The 
lines  show  the  spectra  computed  from  molecular  dynamic$(3)  and  the  dots  show  the  measure¬ 
ments  of  Hanson  and  McTague(7). 
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Wo-h^vt  shown,  for  systems  for  which  the  potential  energy  and  the  appropriate 
connection  to  the  radiation  field  (dipole  moment  or  polarizability)  are  known 
sufficiently  accurately,  that  infrared,  electronic,  and  nonresonance  Raman  spectra  can 
be  computed  from  classical  molecular  dynamics  followed  by  simple  quantum  corrections 
to  the  spectra.  Note  that  no  adjustable  parameters  are  needed  for  any  of  the  spectra 
presented  here.  These  essentially  classical  spectra  can  be  compared  to  experimentally 
measured  spectra  to  check  that  the  underlying  computed  dynamics  are  correct,  and  the 
agreement  illustrated  here  indicates  that  a  basically  classical  view  of  the  atomic  motions 
involved  in  these  spectra  is  a  useful  one,  in  harmony  with  our  well  calibrated  physical 
intuition.  ~ — 
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